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ABSTRACT 

We consider the sensitivity of the circular-orbit adiabatic contraction approximation 
to the baryon condensation rate and the orbital structure of dark matter halos in 
the ACDM paradigm. Using one-dimensional hydrodynamic simulations including the 
dark matter halo mass accretion history and gas cooling, we demonstrate that the adia- 
batic approximation is approximately valid even though halos and disks may assemble 
simultaneously. We further demonstrate the validity of the simple approximation for 
ACDM halos with isotropic velocity distributions using three-dimensional N-body sim- 
ulations. This result is easily understood: an isotropic velocity distribution in a cuspy 
halo requires more circular orbits than radial orbits. Conversely, the approximation 
is poor in the extreme case of a radial orbit halo. It overestimates the response a 
core dark matter halo, where radial orbit fraction is larger. Because no astronomically 
relevant models are dominated by low-angular momentum orbits in the vicinity of the 
disk and the growth time scale is never shorter than a dynamical time, we conclude 
that the adiabatic contraction approximation is useful in modeling the response of 
dark matter halos to the growth of a disk. 
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1 INTRODUCTION 

In hierarchical structure formation, galaxies form in grav- 
itationally collapsing dark matter halos. The dissipative 
baryonic matter cool s and condenses in dark m atter halo 
JWhitefc Reejll978^ . iBlumenthal et al.l Jl98dl described 
the halo response to this condensation assuming a spher- 
ical profile with circul ar orbits and adiabatic disk growth 
JBarnes fc Whitg|l984r) . For an adiabatic change, the angu- 
lar momentum is invariant: J^ oc rM(r) = constant. Then, 
given the distribution of the baryonic disk Md{r) and the 
initial dark matter distribution Mi{r), the final distribution 
of dark matter Mf (r) must satisfy 



Mf{rf)rf ^ M,{r,)r, 



-rrid) 



(1) 



where rud is m ass fraction of the d isk. F urther studies of thi s 
model in clude iRvd en (1988, 1991), and Fl ores et alj Jl993l) . 
Recently I Jesseit. Naab fc Burkert. 1.20021 find that the adi- 
abatic contraction approximation is in agreement with their 
simulations. 

Despite the wide usage of the adiabatic contraction ap- 
proximation, discrepancies between observations and theo- 
retica l predictions motivate a detailed check of its validity. 
First . IBlumenthal et alJ il986l) simply assumes that the disk 
growth time is much longer than the dark matter halo dy- 



namical time. However, individual halos may grow simulta- 
neously with their disks and have different assembly hi stories 
iWechsler et al.ll2002l: IZhao et alJlJOOJlLi et alJlJOOJl . Cos- 
mological simulations of the galaxy formation have shown 
that gas accretion in CDM halos procee ds in t wo distinct 
modes depending on the mass of the halo JKeres et al. 2005| ; 
iBirnboim fc Dekel 1200 j) . In massive halos, gas accretion 
is dominated by "hot mode" in which gas is first heated 
up to the virial temperature of the halo and then cools 
to settle gradually into the halo centre. By contrast, gas 
accretion into small halos is dominated by "cold mode" 
in which cold gas sinks in a dynamical time without be- 
ing shock heated. These scenarios may affect the validity 
of the adiabatic contraction formula. Secondly, dark mat- 
ter is not arranged on cir cular orbits and therefore rM(r ) 



n ot strictly conserved. [ Barnes ! (^83), ISellwoodI il999 
and ISellwood fc McGauehl (I2OO3) report that the approxi- 



mation overestimates th e contraction measu red in simula- 
tions. Moreover, recently iGnedin et al.l i2004l claimed simi- 
lar findings in a cosmological simulation. These authors sug- 
gest that the discrepancy is due to the circular orbit assump- 
tion. 

In this paper, we use idealised numerical experiments to 
investigate the effect of the two assumptions in the adiabatic 
contraction approximation (eq.0 and provide physical intu- 
ition for the numerical trends. In f|H we use one-dimensional 
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Figure 1. The rotation curves of simulated halos with different formation time and gas eooling phases compared with the rotation 
curves with the predictions of adiabatic contraction approximation. The upper-left panel shows the rotation curves for early-formed halo 
{zc = 4) with normal cooling (hot mode). The upper-right panel is for early formed halo (z^ = 4) with artificially high cooling (cold 
mode). The bottom-left panel is for late formed halo {zc = 0.25) with normal cooling, and the bottom-right panel is for late formed halo 
{za = 0.25) with high cooling. In each panel, the diamonds denotes simulation results, the solid lines show the predictions of adiabatic 
contraction approximation, and the dashed lines show the rotation curves of the control model, in which cooling and baryon condensation 
are turned off. 



simulations which incorporate dark matter halo mass accre- 
tion history as well as gas cooling to test the adiabatic disk 
growth assumption in realistic forming halo. We find that 
disk growth time scale is always longer than dynamical time 
of the halo in many cases. We also show that the continued 
dark-matter mass accretion has little affect in the inner halo. 
In ^ we test the circular orbit assumption using high resolu- 
tion N-body simulations with cosmological dark matter halo 
initial conditions. As expected, we find that radial orbits re- 
duce the dark matter halo response to disk growth predicted 
by the simple circular-orbit approach. However, in order to 
maintain an isotropic velocity distribution in cuspy halos, a 
circular population is much larger than radial orbit popula- 
tion. This explains the often-observed consistency between 
simulations of a dark-matter cuspy halo to disk growth. We 
study typical CDM halos in S0|and summarise in ^^ 



2 ADIABATIC CONTRACTION DURING 
CDM HALO FORMATION 

We investigate a dark matter halo response to disk growth 
for several different halo and disk growth time scales using 
one-dimensional hydrodynamic simulations (see Appendix 
A and Lu ct al. 2009). We use 5 x 10* equal- mass shells to 
represent the dark matter distribution and 500 equal-mass 
shells to represent the gas distribution. The evolution of ev- 
ery gas shell is followed until it cools below 10*K; thereafter 
its mass is assigned to an exponential disk. The scale-length 
of the disk, r^, is fixed to ^^^ryjr, where ryjr i s the virial 
radius of the halo at z = (|Mo. Mao, fc WhitelHoflJ^ . The 
disk is assumed to be a rigid exponential disk. 

We examine both an early- {z = 4) and late- (z — 0.25) 
time halo-formation scenarios with M = lO^'^ Mq. Because 
the virial mass of the halo is fixed at 10^^ M©, gas accre- 
tion is primarily in the "hot mode" phase for this halo mass 
JKeres et alJl2005ll . The gas is heated to the virial tempera- 
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Figure 2. The final rotation curve for a circular orbit (long- 
daslied), an isotropic (dotted), and a radial-orbit dominant (dot- 
dashed) distribution. The dark matter halo is c=12 NFW halo 
and disk is fiducial disk. The solid line is for the adiabatic contrac- 
tion approximation prediction. The rotation curve for the adia- 
batic contraction approximation prediction and the circular orbit 
halo case are almost identical. 



ture and then slowly cools to form a disk. In order to examine 
the adiabatic contraction for the "cold mode" gas accretion, 
we artificially increase cooling rate by a factor of 100. Then, 
the accreted gas cools rapidly and joins the disk without 
shock heating. In this case, the disk and its host halo have 
similar growth times. 

For the four simulations, two redshifts and two accretion 
models, we measure the rotation curves at the present time, 
which are shown in Figure^ Since disk growth time is longer 
than dynamical time of the dark matter halo in all four, the 
adiabatic contraction approximation adequately predict the 
response of dark matter halos. Although the host halo con- 
tinues to accrete dark matter, this accretion primarily afi'ects 
the outer halo. In summary, the adiabatic contraction for- 
mula gives an acceptable approximation for these standard 
scenarios. 



3 ADIABATIC CONTRACTION FOR A 

NON-CIRCULAR ORBIT DISTRIBUTION 

We explore disk growth in a dark matter halo represented 
by a Navarro, Frcnk & White (1997, hereafter NFW) pro- 
file, p oc l/[r(r -I- Ts)^] with virial radius r^ir- We as- 
sume a concentration c = r^ir/rs — 12 consistent with 
the rotation curve for a large spiral galaxy and exam- 
ine three cases of different anisotropy to study the ef- 
fects of radial orbits on the circular-orbit adiabatic pre- 
diction: a circular-orbit halo, an isotropic halo, and a ra- 
dial orbit biased halo. The NFW halo with pure circu- 
lar orbits is constructed by assigning each o f the parti- 
cles at radius r a tangential velocity, Vc = \/GM{< r)/r 
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Figure 3. The rotation curve of the pure radial orbit halo. The 
diamonds denote the simulation result, the solid line denotes the 
prediction of adiabatic contraction, and the long dashed line de- 
notes the rotation curve of the reference halo. 



in a random tangential direction. The isotropic and radi- 
ally biased distribution functions are computed using the 
Osipkov-Merritt mo del JBinnev fc TremaindllOSTI : lOsipkovl 
I1979I: iMerritd Il985^ which controls the orbital structure 
by an anisotropy radius ra- For an isotropic halo, Va ~ 
00. For the radially biased case, we choose the minimum 
value for ra that results in positive density. The anisotropy 
profile in this model closel y corresponds to the pro- 
file from a virialized collapse ([Eke. Navarro fc Frenklll998l : 
IColin. KlvDin fc Kravtsol \2Q0(t . The N = 10^-particle 
phase spaces are realized by a Monte Carlo procedure. Our 
rigid exponential disk has mass rrid = 0.04M„ir and the 
disk scale length is r^ = 0.014 r„ir motivated by galax; 
formation in ACDM cosmogony (|Mo. Mao, fc Whitjll99; 
Klvoin, Z hao fc Som crvillc 2 002r . To mimic disk growth, we 
increase the disk mass from zero, keeping the disk scale 
length unchanged. To ensure the validity of the adiabatic ap- 
proximation, the time scale of disk mass growth is 10 times 
longer than the dynamical time of the dark matter halo at 
the disk scale length. The gravitational force on each dark 
matter part icle is calculated using the self-consistent field 
code (SCF^ JClutton-Rrockll97lll97^:lHernuist fc OstrikeH 



ll992l : IWeinberJl999fl . which solves Poisson's equation using 
a set of density-potential bi-orthogonal function expansions. 
Figure |5| shows the total rotation curves after disk con- 
traction for the circular orbit case (long dashed curve), the 
isotropic case (dotted curve) and the radial orbit biased case 
(dot-dashed curve). For comparison, we also show the pre- 
diction of the adiabatic contraction approximation (solid 
curve). The adiabatic contraction approximation is nearly 
exact for circular orbits, but overestimates the contraction 
when the eccentric orbit contribution increases. The radial 
orbit biased halo model is not significantly biased toward 
radial orbit around disk (ra = O.lrvir — 7rd)- Consequently 
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Figure 4. Rotation curves from simulations cooperated with the adiabatic contraction approximation for isotropic NFW and core halos. 
The dotted line is the rotation curve from disk only, the solid line is the total (disk + dark matter halo) rotation curve from the simulation 
results, and the long dashed line is the total rotation curve from the a diabatic contraction a pproximation. We also plot the total rotation 
curve predicted by the modified adiabatic contraction approximation of lGnedin et al.N2004) with A = 0.85 and w = 0.8 (dot-dashed line). 
The disk has the fiducial parameters: r^ = 0.014r„i,. and m^ = 0.04m„ir- The overestimation by the adiabatic contraction approximation 
increases as halo concentration decreases. The modified adiabatic contraction approximation provides a much better fit to the simulation 
results. 



the effect of radial orbit is not dramatic in Figure |5| How- 
ever it is clear that the radial orbit reduces the dark matter 
halo response. 

Since it is difffcult to simulate pure radial orbit halos 
in three-dimensional N-body simulation, we use implement 
one-dimensional simulation with a pure radial orbit halo. 
The orbits are purely radial and the density profile of the 
simulated halo is proportional to r~^ in central region. Fig- 
ure 13 shows the rotation curve of the simulated halo com- 
pared with the prediction of adiabatic contraction formula. 
The circular-orbit prediction overestimates the contraction 
by factor of two. 



4 ADIABATIC CONTRACTION FOR 
TYPICAL CDM HALOS 

We consider a range of halo parameters and disk masses 
to explore the general applicability of the circular- 



orbit adiabatic approximation. Recent cosmological simu 
lations show that velocities in the inner region o f dark 
matter halos is isotropic llEke. Navarro fc FrenkI 






1996 



2001 



[Colin. Klvpin fc Kravtsov"2000': 'Fukus hige fc Makinol 
[picmand, Mooro fc Stadcl 2004). Therefore we explore 
isotropic halos. Although an NFW halo model is cur- 
rently accept ed in CDM cosmology, some recent theoret- 
ical models llMo fc Mad l200l l2004t lOh fc Bensonl l2003l: 
'Weinberg fc Kat3 |2002) and observations JDe Blok et alJ 
2001) suggest that dark halos may have cores. Simulations 
for the core halo are also carried for the three cases, with 
c = Tyir/rcore = 15, 12, and 7.5, respectively. In addition 
to the fiducial disk, mrid — 0.04Mj,ir and r^ = 0.0147?„ir, 
we consider a low-mass [nid ~ 0.02M„ir) and a high-mass 
{rrtd = O.lMvir) disk. 

Figure 2] compares the post-formation rotation curves 
for NFW halos and core halos with c = 15 and c = 7.5 with 
the adiabatic contraction predictions. The formula overes- 
timates the rotation velocity for these astronomically mo- 
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tivated halo models although the discrepancy is modest. 
Figure 13 quantifies the relative differences, r) = \Vad — 
Vsim\/Vaim, as a fuuctiou of halo concentration c, where 
Vsim and Vad are the circular velocities at r = 2.2Rd ob- 
tained from simulation and from the adiabatic contraction 
(eq. , respectively. The discrepancy is 4% for NFW model 
with c = 15 and increases to 8% for c = 7.5. The discrepancy 
increases with disk mass but the dependence is weak. The 
c = 7.5 case is a low value for galaxy halos in the current 
CDM model, therefore, these results show that the adiabatic 
contraction approximation remains good for isotropic NFW 
halos. For core halos, the discrepancy is 14% for c = 15 and 
increases to 23% for c = 7.5 with weak dependence on md- 



iGnedin et al.l ll2004f) improved the agreement between 
the adiabatic contraction approximation and the simulation 
result using a modified version of the original adiabatic con- 
traction approximation including the effect of non-circular 
orbits. In this model, the adiabatic invariant M{r)r is re- 
placed by M(f)r, where r = Ar™ . The authors suggested 
A ~ 0.85 ± 0.05 and w f» 0.8 ± 0.02 based on cosmological 
dark matter halo simulations. As comparison, we show in 
Figure |1] the rotation curves obtained from this modified 
adiabatic contraction approximation, together with those 
obtained from the simulation and from the original adia- 
batic contraction approximation. As one can see, the mod- 
ified adiabatic contraction approximation agrees with the 
simulation results for NFW halos. For halos with a con- 
stant density core, there is clear discrepancy between the 
model and simulation, although the modified model is better 
than the original ad iabatic contraction model. Although the 
iGnedin et alJ l|2004h model takes the effect of non-circular 
orbits into account, their fitting formula, r = Ar^ , is not 
based on the orbital structure. It simply reduces the halo 
contraction using fitting formula with empirically su ggested 
fitting parameters. Therefore the lGnedin et alJ l|200^ model 
estimation for core halo cases with suggested fitting param- 
eters does not work as well as one for NFW halo cases. For 
core halo, different fitting parameters and further modifica- 
tion in fitting formula are required. 

Figure 1^ shows that the adiabatic contraction approx- 
imation is better for high-concentration halos. This trend 
is explained by the distribution of orbits. We describe the 
shape of an orbit by the ratio of the angular momentum to 
the angular momentum of a circular orbit at a fixed energy 
K — J/Jmax{E). The orbit with k = (k = 1) is radial (cir- 
cular). Figure |S| describes the ensemble average of this ratio 
for isotropic halo in radial bins (K)(r). At a fixed energy, the 
mean value of k is 2/3 for comparison. 

Figure |S| shows that an isotropic NFW halo has more 
low-eccentricity orbits at fixed radius than an isotropic core 
halo. This can be understood as follows. The density at a 
given radius in a halo is contributed by particles on differ- 
ent orbits. In the inner region of a halo, orbits with lower 
energies are more circular, while those with higher energies 
are more radial. Assuming isotropic velocity dispersion, one 
can show that the energy distribution is flatter in a core 
halo than in a cuspy halo. Consequently, for isotropic veloc- 
ity distribution a core halo requires more high-eccentricity 
orbits than an NFW halo. 
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Figure 5. Relative difference (rj = °'''° 
rotation speed from the simulation and from the adiabatic con- 
tract approximation at 2.2Rd for c=15, 12, and 7.5 NFW and 
core halos. Three disks, m^ =0.1, 0.04, and 0.01 Myi^, are con- 
sidered for each halo. The relative difference depends strongly on 
the halo structure and the disk mass dependence is negligible. 
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Figure 6. The mean relative angular momentum per orbit, k = 
J/Jmax{E), is averaged in radial bins for both the NFW and 
core halo models and three concentrations. The top three lines 
represent NFW halo models and bottom three lines represent 
core halo models. The solid lines show c=15 NFW and core halo 
models, the dotted lines shows c=12 halo models, and dot-dashed 
lines shows c=7.5 halo models. The fraction of eccentric orbits 
increases dramatically in core halos. 



5 SUMMARY 

We study th e accuracy of the circu lar-orbit adiabatic ap- 
proximation JBlumenthal et alJl986 l) in predicting halo con- 
traction due to disk formation and provide a physical expla- 
nation for the deserved trends. We consider: (1) variation 
in the accretion time scale; (2) variation in the accreted 
disk mass; (3) variation in the central concentration of both 
cuspy and core halo models; and (4) variation in the velocity 
isotropy. The circular orbit adiabatic contraction approxi- 
mation is acceptable over a wide range of astronomically 
interesting parameters. The relative change in the rotation 
curve value between the simulation and circular orbit ap- 
proximation at 2.2 disk scale lengths, rj, is less than 23% 
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for our entire range of realistic parameters. We find tfiat 
the disk growth is still slower than dark matter halo dy- 
namical time in the vicinity of the disk and therefore the 
adiabatic approximation is maintained. The value of rj de- 
pends only weakly on the fraction of accreted baryon mass, 
and therefore, the dependence on halo concentration cannot 
be explained by disk dominance in less concentrated halos. 
However, 77 is strongly correlated with the fraction of ec- 
centric orbits in the distribution. The steeper the cusp, the 
larger fraction of more circular orbits are required at fixed 
radius, and this supports the circular orbit approximation. 
Although the adiabatic concentration approximation over- 
estimates the response of dark matter halos, as long as the 
dark matter halo has central cusp and isotropic velocity dis- 
tribution the overestimation is negligible. 

Our results have important implications for the forma- 
tion of disk galaxies in the CDM scenario. In comparing 
theory with the observed TuUy-Fisher relation, one usu- 
ally uses the peak rotation velocity of galaxy disks to rep- 
resen t the observed rotatio n velocities iIMo. Mao, fc White! 
Il998t) . However, as shown in lMo fc Mad ll200(f) . if dark mat- 
ter halos respond to the disk growth according to the adi- 
abatic contraction model, current CDM model predicts a 
TuUy-Fisher relation that has a much too low zero-point 
(i.e. galaxies are too faint for a given peak rotation ve- 
locity). In order to match the observed TuUy-Fisher zero- 
point, one has to assume that disk growth d oes not cause 
any c ontraction in dark matter halos at all (jCroton et alJ 
I2OO6I . e.g.). This assumption is not supported by our re- 
sults, which show that adiabatic contraction approxima- 
tion works reasonably well for CDM halos over a wide 
range of situations. However, it should be pointed out 
that there are other astrophysical p r ocesses, such as merg- 



ers (Dekel. Devor. fc Hetzronil 12003: iBovlan -Kolc hin fc Ma 



2004), dynamical heating by substructure s (El-Z ant et al 



2004), halo pre-processing ||Mo_fc_Map '2004), and resonance 
dynamics llWeinberg fc Katal2002f) . that may modify halo 
structures but are not included in the models considered 
here. Unfortunately, the importance of these processes are 
not well understood at the present. 
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APPENDIX A: ONE-DIMENSIONAL 
HYDRODYNAMIC SIMULATIONS 

We use the Lagrangia n based one-dim ensional hydrody- 
namic code described in lLu et alj i200dl to simulate forma- 
tion of a halo from an initial perturbation at a high redshift 
Zj. We assumes the halo mass accretion histories proposed 
bv lWechsler et alJ ^2002^ 



M{a) = Mo exp 



--•(?-') 



(Al) 



where a is the expansion scale factor, Oc is the scale factor 
corresponding to the formation time of the halo, and Mo is 
mass of the halo at the observation time ao. In this func- 
tion, ttc is the only free parameter to characterize the shape 
of a m ass accretion history. Reader may refer to ILu et alJ 
i200tj) for detailed description on making the initial condi- 
tions given the mass accretion history. We choose Zi = 200, 
and the initial temperature of gas shells is set to be the CMB 
temperature at this epoch. 

The simulation has both dark matter shells and gas 
shells. The gas initially follows the distribution of the dark 
matter but evolves differently from the dark matter due 
to hydrodynamics. We use the Lagrangian finite-difference 
scheme to follow the evolution of the shells. The numerical 
treatm ent is same as what is described inlThoul & Weinberg 
il995h . To avoid numerical instability due to dark matter 
shell-crossing, the mass of each dark matter shell is chosen 
to be much smaller than that of a gas shell. The baryon 
fraction is fixed at /(, = 0.17. The chemical abundance is 
assumed to be primord i al. Th e radiative cooling function 
proposed bv lKatz et alJ il996l) is implemented in the simu- 
lations. 

When a gas shell cools to a temperature, below 10*K, 
the gas in the shell is considered to be cold. Since we do not 
include any cooling processes below this temperature, the 
cold gas is assumed to retain a temperature of IQ^K until it 
flows into the center of the halo. At this point, the gas joins 
a central exponential disk with a scalelength r^ — 2i^rvir, 
where r-vir is the virial radius of the halo at z = 0. The cold 
gas disk is assumed to be a rigid object, and its gravity is 
included in the subsequent evolution of other mass shells. 
At any given time, the gravitational acceleration of a shell 
at radius r^ is given by 

„2n GM(< n)ri 

g. = i/of^AH - ^^2 ^ ^2p/2 , (A2) 

where Hq is the Hubble's constant at the present time, f^A is 
the density parameter of the cosmological constant, M(< r;) 
is the total mass (dark matter, gas and exponential disk) 
enclosed by r;, and a is a softening length taken to be 0.0005 
of the virial radius of the halo at the present time. This scale 
is much smaller than any scale of interest. In the simulations, 
the density parameter of the non-relativistic matter and of 
the cosmological constant are JIm = 0.3 and JIa = 0.7. and 
the Hubble's constant is ffo = 100km s~^ Mpc~^. 



